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7-rays from annihilating dark matter in galaxy clusters: stacking vs 
single source analysis 
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ABSTRACT 

Clusters of galaxies are potentially important targets for indirect searches for dark matter 
annihilation. Here we reassess the detection prospects for annihilation in massive halos, based 
on a statistical investigation of 1743 clusters in the new Meta-Catalog of X-ray Clusters. We 
derive a new limit for the extra-galactic dark matter annihilation background of at least 20% 
of that originating from the Galaxy for an integration angle of 0.1°. The number of clusters 
scales as a power law with their brightness (boosted by dark-matter substructures), suggest- 
ing that stacking may provide a significant improvement over a single target analysis. The 
mean angle containing 80% of the dark-matter signal for the sample (assuming an NFW DM 
profile) is ~ 0.15° (excluding the contribution from the PSF of any instrument), indicating 
that instruments with this angular resolution or better would be optimal for a cluster annihi- 
lation search based on stacking. A detailed study based on the Fermi-LAT performance and 
position-dependent background, suggests that stacking may result in a factor ^2 improve- 
ment in sensitivity, depending on the source selection criteria. Based on the expected perfor- 
mance of CTA, we find no improvement with stacking, due to the requirement for pointed 
observations. We note that several potentially important targets: Opiuchius, A 2199, A 3627 
(Norma) and CIZA J1324.7— 5736 may be disfavoured due to a poor contrast with respect to 
the Galactic dark-matter signal. The use of the homogenised MCXC meta-catalogue provides 
a robust ranking of the targets, although the absolute value of their signal depends on the exact 
dark matter substructure content. For conservative assumptions, we find that galaxy clusters 
(with or without stacking) can probe (erf) down to 10^^^ — 10^^"* cm"^ s^^ for dark mat- 
ter masses in the range 10 GeV - 100 GeV. For more favourable substructure configurations, 
{av) ^ 10"^® cm'^ s~^ may be reached. 
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1 INTRODUCTION 

The annihilation of dark matter (DM) particles into 7— rays has 
been flagged as one of the most promising channels for in- 
direct detection. Regions of high DM density are of particu- 
lar interest, making th e Galactic centre the most obvious target 
JSilk & Bloemenll987h . However, the Galactic centre is plagued by 
a large astrophysical 7— ray background at all angular scales that 
make s any DM signal difficult to identify (e.g., [Aharonian et al.l 
1 20041) . In that respect, dwarf spheroidal galaxies (dSphs) have 
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the advantage to be essentially background-free, relatively close 
by and with DM density profiles that can be constrained from 
their internal kinematics . This has made th e m popular candidates 



for i n direct detectio n ( Evans et al. 20041: Bergstrom & Hooped 
20061: IStrigari et al.l l2007t iPieri et al I I2OO9I: |Ab 



Kuhlen '2OIO': 'Sanchez-Conde et alj 120111: IWalkeretalJ 120111 : 
iCharbonnier et al. 201 1 : Acker mann et al.ll201 ill 

Somewhat less explored to date, clusters of galaxies are the 
largest gravitationally bound structures in the universe, the large 
DM content of which ma kes them potentially inte resting tar- 
gets for indirect detection JColafrancesco et al.l 120061) . Although 
strong constraints have already been derived from X-ray and 
gravitational tensing stud i es on the DM distribution in clusters 
jPointecouteau et alj20o3 : IVikhlimn et al.l2006l : lBuote et alj2007l : 
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IShan et alj|2olol : IPastor Mira et alJlioTTI: lEttori et al.ll201 ih . con- 
straining the inner DM distribution is still a challenging task. 
Even strong lensing, which likely is the best suited way to pin 
down the DM distribution at the cluster centre, fails to assem- 
ble convincing cons traints (see for instance the different con- 
clusions reached by Limousin et alj uOOTt iNewman et alj uOllI : 



iMorandi & LimousinI 



20121) . Estimates of the DM profile and cal- 



culations of the 7— ray flux fro m clusters are based on X-ray ob- 
servations, fro m which NFW ([Navarro. Frenk & Whitel 119971) or 
Einasto (e.g.. iMerritt et al.l l2006h profiles are assumed. For in- 
stanc e, based on the HIFLUGCS catalogue contain ing 106 ob- 
jects JReiprich & B6hringenl2002t IChen et al.ll2007h . several au- 
thors have identified the potentially most lumino us objects in 
DM emission, such as Fornax, Coma or Perseus jjeltema et al.l 
12009 ; Pinzke et al. 2()ir). The non-detection of these favoured 
targets by Fermi-LAT and H.E.S.S. has resulted in constraints 
on the DM annihilation cross-section 1 Ackermann et al.l 20 id : 
lYuan et all l20ld; lAndo & Nagail 1201 2l : lAbramowski et al.ll2012h . 



See, however. 



LU; lAndo & JNaga 
Han et all ( l2oIa) ^ 



for a possible evidence of an ex- 
tended emiss i on. A lternatively to these ' observational' approaches, 
ICuesta et alj 1 1201 11) have performed synthetic Fermi observations 
from the CLUES constrained cosmological N-body simulation of 
the local universe and flagged Virgo and Coma , along with DM 
filaments as interesting targets. lOao et al.l ( 120 12h is another exam- 
ple of high-resolution N-body simulations used to estimate the DM 
profile/content and signal of selected targets. 



In this study we make use of t he recently published Meta- 
Catalogue of X-ray Clusters, MCXC JPiffaretti et ai1l201lh . which 
contains 1743 clusters of galaxies. The size of the catalogue, with 
~ 17 times more objects than the HIFLUGCS catalogue, makes 
it possible to investigate some statistical aspects of DM indirect 
detecti on in galaxy cluster s. This paper is part of a series: a first 
paper dCombet et al.l 120121) highlighted the improvement brought 
by a stacking analysis over a single source analysis for the DM de- 
cay case. The current paper focuses on the DM annihilation case: 
we provide a quantitative analysis of the best observing strategy to 
use for the Fermi-LAT and CTA observatories, we discuss the po- 
tential benefit of a stacking strategy with respect to single source 
observation, and we also present the number of objects to look at 
to optimise detectability. The last paper of the series addresses the 
possibility of using the sta cking analysis to dis entangle CR-induced 
from DM-induced signal (iMaurin et al.ll2012r) . 



2 THE MODEL AND ITS INGREDIENTS 

The 7-ray flux <l>.y from dark matter annihilations 
(cm^^ s^^ sr^^ GeV^^) received on Earth in a solid angle 
Afi, is given b)[j 



IK 



1 (o-annf) 

47r (5m I 






(1) 



where 5 = 2 for a self-conjugate particle and 4 otherwise, m-^ is the 
particle mass, (aann^') is the velocity-averaged annihilation cross 
section, dN-,/dE~, is the energy spectrum of annihilation products. 



2.1 Spectrum and astrophysical factor J 

The differential annihilation spectrum, dN.,/dE^, requires a spe- 
cific DM particle model. It is the sum of a prompt contribution and a 
contribution from inverse Compton scattered (I CS) secondary elec - 
trons and positrons with the CMB (see, e.g., iHuang et alj|2012h . 
For the sake of simplicity and to keep the analysis as DM par- 
ticle model-independent as possible, we disregard the 'delayed' 
ICS contribution. The latter has a si milar spatial distribution to 
that of the prompt dHuang et alj|2012i) . so that the factorisation of 
the spatial and energy-dependent term in Eq. ([T} holds. Actually, 
depending on the annihilation channel, the ICS contribution can 
dominate over the prompt one. Considering only the prompt con- 
tribution as we do here provides a conservative and robust lower 
limit on detectability. In this paper, we further restrict ourselves 
to the bb annihilation chan nel, taken from Eq. (6) and Table XXII 
l20Ilh . We note that the spectral parameters 
201 lb are provided for WIMP masses in the 



Cembranos et al.l 



Cembranos et al.l 



m 
in 

range of 50 GeV to 8 TeV. Here we assume the spectral parame- 
ters for masses below 50 GeV are given by the parameters for a 50 
GeV mass, and similarly above 8 TeV. The results are not strongly 
affected (less than a factor 1 .5 in the sensitivity limits) by the choice 
of the 7-ray annihilation channel (apart from the tt channel). 

The ' J-factor' represents the astrophysical contribution to the 
signal and corresponds to the integral of the squared dark matter 
density, p^il, Jl), over line of sight I and solid angle Af2, 



J(At7) = / / p''{l,n)dldn. 
'An . 



(2) 



We have AJl = 27r ■ (1 — cos(Qint)), and aint is referred 
to as the 'integration angle' in the following. All J-factors pre- 
sented below, including substructures (in the Galaxy or in galaxy 
clusters), are calculated from the public code CLUMPY v20 11.09 
dCharbonnier et al.ll2012l) . 



The paper is organised as follows: in Section |2] we briefly 
present the key quantities for the signal calculation (J-factor, DM 
halo profiles). In Section|3] the MCXC catalogue is introduced, and 
the cluster signal distribution presented, along with the resulting 
skymap. The contrast with the Galactic DM annihilation signal and 
the astrophysical background, and the consequences for the ranking 
of the best targets are also discussed. The stacking approach and re- 
sults are presented in Section |4] In particular, the boost of the DM 
signal from DM substructures (in the galaxy clusters) and its effect 
on the stacking is detailed. The sensitivity to a DM signal taking 
into account realistic instrumental responses is then evaluated for 
Fermi-LAT and CTA instruments. We conclude in Section |5] (Ap- 
pendix |A]provides parametric formulae to evaluate the signal from 
a cluster for any integration angle. Appendix |B] provides a quick 
comparison to values of J found in other works). 



2.2 The smooth DM halo and substructures 

For the DM halo smoot h profile, we use an NFW 
dNavarro. Frenk & Whitel 1 997h 



p{r) 



fe)('-t)- 



(3) 



^ We remind that the spatial term J in Eq. {T) couples to the energy- 
dependent temi dN-y/dEj for objects at cosmological dista nces , because 
7-rays are absorbed along the line of sight (e.g.,[Cirelli et al. 2 010l) . Th e red- 
shift d istribution of the MCXC catalogue of galaxy c lusters (P iffaretti et aj 
I2OI ih peaks at ^ ~ 0.1 (see their Fig. 1): following ICombet et al. t2012r 
we neglect the absorption for the MCXC galaxy clusters. 
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where r^ is the scale radius and ps is the normaHsatior|j. We note 
that Einasto profiles give slightly more 'signal' than NFW halos, 
making our conclusions on detectability conservative. 

Cold DM N-body simula tions show a high leve l of dumpiness 
in the DM distribution (e.g., iDiemand et al.ll2007l : ISpringel etal] 
120081) . These substructures boost the signal i n the outer parts o f the 
DM halos. In agreement with the analysis of ICao et al.l ((2013), we 
find that the boost in galaxy clusters is larger than the boost ob- 
tained for less massive objects suc h as dSphs. For the latter, boost 
are < 2 JCharbonnier et al.ll201 lb . whereas we obtain an overall 
boost of ~ 10 — 20 for galaxy clusters based on conservative as- 
sumptions for the substructure parameters (the impact of these pa- 
rameters in discussed in Section U!4t . The reason is twofold: first, 
dSphs are less massive so that the mass range of substructures is 
smaller (the minimal mass is assumed to be the same regardless of 
the object), hence the number of objects, and their overall contribu- 
tion; second, the effective angular size on the sky is larger for dSphs 
so that current instrumen ts integrating ou t to 0.5° integrate less sub- 
structure signal (see also lOao et al.l2012h . These boost are obtained 
from the following configuration — used throughout the paper with 
the exception of Section [J!4] — for the mass and spatial distribution 
of the substructures : i) dNsubs/ dM oc A/~^'^ with a mass frac- 
tion / = 10% in substructures JSpringel et al.luOOSi) . a minimal 
and ma ximal mass of 10 ~^ Mq and lO^'^Mdustcr respectively, 
and the iBullock et alj ( 1200 If) concentration (down to the minimal 
mass); ii) the substructure spatial distribution dNsubs/dV follows 
the host halo smooth profile. For this configuration, we checked 
that the boost is only mildly dependent on this mass by varying 
the mass from lO^'' to 1 Mq. We note that the minimal mass for 
sub-halos can be a s small as 10~^° Mg depending on the particle 
physics model (see lProfumo et al.ll200q , and references therein). 

A complete study of the boost should consider different pro- 
files, different parametrisations for the mass-concentration relation- 
ship, etc. This will be fully addressed in a future work. However, 
given the impact it can have on the {(Tanni') limit (or detectability 
for current instruments), a short discussion and general trends are 
given in Section [J!4l 



3 J FACTORS FOR THE MCXC SAMPLE 

The MCXC dPiffaretti et ai]|201ll) contains 1743 clusters of galax- 
ies detected in X-rays, and assembled from publicly available cat- 
alogues mainly based on the ROSAT All Sky Survey or ROSAT 
serendipitous catalogues. Most observational constraints and pre- 
dictions are expressed in terms of A = 500 or A = 200. For in- 
stance, the mass of a halo. Ma can be defined within a radius Ra 
within which the average density reached A times the critical den- 
sity of the Universe (at a given redshift). The MCXC provides ho- 
mogenised quantities for each clusters computed within A = 500, 
e.g., the standardised [0. 1-2.4] keV X-ray luminosity L500, the total 
mass Afsoo, the radius -Rsoo- 

To fully describe the NFW profile parameters (see Eq. |3] for 
each galaxy cluster of the MCXC catalogue, we used the provided 

^ A decreasing inner slope with the halo radius r (Einasto profiles) tends 
to be f avoured by recent h i gh-resolution N-body simulations ( NavaiTO e t al.| 



|2004 iMerritt et al.l I2OO6I: ISpringel et al."2008 yMaitizzi et alj 120121) and 
also by galaxy observations (Cheminetal. 2 111) . Simulations including 
baryons and feedback processes are important to address fu rther the ques- 
tion o f the (dark) matter profile in the innermost region (e.g.. lMartizzi et alj 
l20il . 



M500 together with a mass-concentration relationship (i.e., ca is 
fully determined by the cluster mass A/a). This relation is ob- 
serva t ionally constrained at the cluster sc ale (IPointecouteau et al.l 
l2005l ; lBuote et al.ll2007l ; lEttori et al.ll2010l) . It has also been shown 
to depend on the epoch of ha lo formation by numerical sim- 
ulations of structure formation C Bullock et al.l BOOU iDolag et al.l 
I2OO4I ; iDuffv et alj [20o3 ; iKlvpinetal. 20nir Although the data 
present a large dispersion , a systematic offset remains unexplained 
Duffv et al.i2"ooi . l2010h . In this study, we assume the lDuffv et ahl 



20081) mass-concentration relation. 



For an NFW profile Vs — Ra/ca and the scale density ps is 
obtained from the mass measurement. The J factor for all clusters 
are then calculated from Eqs. (|2) and (|3} with CLUMPY. 



3.1 Brightest targets 

Figure [T] provides a synthetic view of the J-factor for each galaxy 
cluster of the MCXC catalogue as a function of their angle (f} 
away from the Galactic centre. The integration angle is taken to 
be Oint = 0.1° (left panel) and aint ~ 0.5° (right panel), the 
typical range of value for the energy-dependent angular resolu- 
tion of current 7-ray instruments such as Fermi-LAT (in the high- 
energy range above ~10 GeV) and H.E.S.S. Table|T]gathers results 
for the twenty brightest clusters in the MCXC. From this table, 
we simply note that J-factors are c ompetitive with those obtained 
for dSphs (e.g.. IWalker et al.ll201 ih . confirming that galaxy clus- 
ters are valid targets for dark matter annihila tion searches (see also 
ISanchez-Conde et al.ll2oTTI ; lOao et al.ll2012h . The panels of Fig.|2] 
show a skymap version of Fig. [T] The top left panel shows the J 
factor induced by DM annihilation in the Galactic halo cumulated 
with all MCXC objects. The top right panel shows the J factor 
skymap for all MCXC galaxy clusters only. The bottom panel lo- 
cates the twenty most promising targets labelled by distance, ab- 
solute J-factor value and contrast with respect to the DM Galactic 
signal. 

Several clusters including Virgo, Coma, Fornax, NGC 5813 
and Ophiuchus have already been credited to be interest- 
ing sources in numerous stu dies give n their mass e s and 

distances dColafrancesco et al.l l200q; Ijeltema et al.l 200S; 



Sanchez-Conde et al.l I2OI ll ; IPinzke et al.l I2OI ll ; ICao et alj 12012 : 
Han et all l2012h . Other objects, such as, 3C129 and AWM7 



were only highligh t ed from the HIFLUGCS catalogue a nalysis 
djeltema et alli2009l ; IPinzke et alj|201 ll ; iHuang et aljhod) . With 
ten times more objects, the MCXC gives a more exhaustive list of 
potential targets including, e.g., J0123.6-I-3315 and J1324.7-5736 
(see Table[T]and Fig.|2]l. 

Some differences exist with previous calculations (see 
App. |BJ. These can be partly attributed to a different prescrip- 
tion for the substructures. However, another important difference 
comes from the fact that almost all previous studies are based 
on the Afsoo va lues obtained from the HI FLUGCS catalogue 
dReiprich & Bohr inger 2002; C hen et al.l2007h . In particular, some 
of 'brightest' objects found (e.g.. Coma, Fornax, AMW7) have 
larger masses than thos e provided in the MCX C catalogue. As dis- 
cussed in the App. A of lPiffaretti et alJd20lil) , the MCXC relies on 
a more accurate model for the gas distribution, and many compar- 
isons to numerical simulation s indicate that any systematic uncer- 
tainties are now < 15 - 20% dPiffaretti & Valdaminill2008l) . 
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Figure 1. Computed J-factors for the MCXC sources (the 10 highest-contrast clusters are hi ghlighted, the remaining a re shown with a '+' symbol) vs Galactic 
DM background (total is the sum of smooth, sub-halos, and cross-product — see details in ICharbonnier et al.ll2012h . The yellow filled square symbols are 
evaluated from the cumulative of the cluster signal in different (p bins: this can be interpreted as a lower limit for the extra-galactic DM annihilation signal. 
Left panel: integration angle a^^t = 0.1°. Right panel: ctint = 0.5°. 



Table 1. Twenty brightest galaxy clusters from the MCXC and their contrast J/ Jcal f™' CKin 
position of the cluster (angle ij> away from the Galactic centre, see Fig.[T). 



0.1°. The DM Galactic background is evaluated at the 



Name 



Index / 

MCXC (deg) 



(deg) (deg) 



d 
(Mpc) 



lof 



[rank] ' 



10 (H^) ^80% logio [JKnt)/(M| kpc-5M ^(£ 

- " (deg) (0.1°) (0.5°) (080%) (0.1°) (0.5°) 



t) 

JGal("int) 



Virgo 

A426 

A3526* 

NGC 4636 

A3627*-^ 

Coma 

NGC5813* 

Ophiuchus* '" 

NGC5044* 

AWM7 

A1060 

Fornax 

A 1367 

J 1324.7-5736* 

A0262 

3C129'' 

A2199* 

NGC1550 

A3S71* 

2A0335 



884 


283.8 


74.4 


86.3 


15.4 


14.3 


258 


150.6 


-13.3 


148.0 


75 


15.1 


915 


302.4 


21.6 


60.1 


48.1 


14.5 


906 


297.7 


65.5 


78.9 


13.2 


13.3 


1231 


325.3 


-7.1 


35.4 


66 


14.6 


943 


57.2 


88.0 


88.9 


96.2 


14.9 


1147 


359.2 


49.8 


49.8 


21.3 


13.6 


1304 


0.6 


9.3 


9.3 


116. 


15.0 


978 


311.2 


46.1 


62.8 


36.9 


14.0 


224 


146.3 


-15.6 


143.3 


72.1 


14.5 


689 


269.6 


26.5 


90.4 


53.1 


14.2 


285 


236.7 


-53.6 


109.0 


21.7 


13.5 


792 


235.1 


73.0 


99.6 


89.3 


14.6 


990 


307.4 


5.0 


52.8 


79.5 


14.5 


158 


136.6 


-25.1 


131.1 


68.4 


14.3 


350 


160.5 


0.3 


160.5 


91.7 


14.5 


1249 


62.9 


43.7 


70.8 


12.4 


14.7 


324 


191.0 


-31.8 


146.5 


55.2 


14.1 


1048 


316.3 


28.6 


50.6 


159.7 


14.9 


286 


176.3 


-35.5 


144.8 


142.5 


14.8 



3.3 
1.2 
1.2 
1.7 
0.9 
0.8 
1.4 
0.7 
1.0 
0.8 
0.9 
1.2 
0.7 
0.7 
0.7 
0.7 
0.6 
0.8 
0.5 
0.5 



11.1 
10.8 
10.7 
10.6 
10.6 
10.6 
10.6 
10.6 
10.5 
10.5 
10.5 
10.5 
10.5 
10.5 
10.5 
10.4 
10.4 
10.4 
10.4 
10.4 



11.8 
11.5 
11.4 
11.4 
11.3 
11.3 
11.3 
11.2 
11.2 
11.2 
11.2 
11.2 
11.1 
11.1 
11.1 
11.1 
11.0 
11.1 
11.0 
11.0 



12.6 
11.8 
11.7 
11.8 
11.5 
11.5 
11.6 
11.4 
11.5 
11.3 
11.4 
11.5 
11.2 
11.3 
11.3 
11.2 
11.1 
11.2 
11.0 
11.0 



20.7 [2] 
21.2 [1] 

4.7 [30] 

6.8 [13] 

1.4 [-] 
7.7 [10] 
2.7 [-] 
0.1 [-] 

3.6 [-] 
11.5 [3] 

6.7 [14] 
8.6 [7] 

6.9 [12] 
2.3 [-] 
9.0 [6] 

10.2 [4] 

3.5 [-] 
9.2 [5] 
1.9 [-] 

8.6 [8] 



4.9 [1] 

4.5 [2] 
0.9 [17] 
1.4 [9] 
0.3 [-] 

1.6 [8] 
0.6 [39] 
0.02[-] 
0.7 [-] 
2.3 [3] 
1.3 [11] 

1.6 [7] 

1.3 [12] 
0.4[-] 

1.7 [6] 

1.8 [4] 
0.6 [38] 
1.7 [5] 
0.3 [-] 

1.4 [10] 



' Whenever the rank is larger than 50, we use [-]. 

* Weakly contrasted clusters are probably not the best targets. 

* Clusters close to the Galactic plane are not favoured targets. 



3.2 Galactic and extra-galactic DIVI background 

Galactic DM provides a 'diffuse' DIVI emission Jcai that can drown 
the point-like emissions we are looking for. The value of the local 
DM density is still loos ely constrained in [ 0.2 — 0.4] GeV cm~'^ 
by several techniques (Sofue et al.l |2009|: ICatena & Ulliol |2010| ; 
ISalucci et allbOld - lOarbari et al.ll201 ll: llocco et al.ll201 ih . We as- 
sume here pq — 0.3 GeV cm~^. The value for JGa.i(4> ^ 20°) is 
also very sensitive to the Galactic sub-halo distribution. The Galac- 
tic signal is thus uncertain by a factor of a few. We calculate in 
Table[T]the contrast, i.e., the ratio between the cluster signal to the 



DM Galactic signal. As shown in Figs. \T\ and [2] the DM Galac- 
tic signal has a shallow latitudinal dependence except towards the 
Galactic centre (6 > 5°) where the signal is maximal. Several of 
the brightest sources are close to the galactic centre, namely Ophi- 
uchus, A3627(Norma), and J1324.7-5736. Although they exhibit a 
large J-factor, their contrast is low, and they are not favoured. In- 
deed, the contrast indicates when a point-like observation strategy 
becomes less promising than a strategy based, e.g., on the detection 
of a gradie nt for smooth Galactic ha lo towards the Galactic centre 
(as done in lAbramowski et al.ll201 ih . Away from the Galactic cen- 
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Figure 2. Top panels: J-factor skymap for aint = 0.1° for Galactic + MCXC sources (left) and MCXC sources only (right). Bottom panel: positions of the 
20 closest (red circles), biightest (black points), and highest J/ Joal (blue circles) from the MCXC. 



tre, we have Jcai oc aint- This is illustrated by the left and right 
panels of Fig.[T] where the value of Jcai is multiplied by 25 mov- 
ing from aint ~ 0.1° to aint ~ 0.5°. However, the corresponding 
signal from each cluster is only marginally increased, meaning that 
the contrast is worsened for large integration angles. 



The diffuse extra-Galactic DM signal constitutes another 
background, the level of which has been estima ted from N-body 
simulations (see, e.g.. Fig. 4 of IPieri et al.ll201 ll) . It is not consid- 
ered here. However, by averaging in each <j) bin the signal from all 
clusters and correcting for the solid angle element, we derive a first 
'data-driven' estimate of this extra-galactic contribution, and we 



find J< 



extra — gal 



Jcai/S (yellow filled squares in Fig.[TJ. Larger 



samples of galaxy clusters are required to refine this figure. 



The five brightest sources in Table [T] are located far from the 
Galactic centre and plane, and therefore have the 'best' contrast 
w.r.t. diffuse DM and astrophysical emissions (located mostly in the 
disk). These sources are also amongst the closest targets and have 
the largest angular size. As we will show in the next sections, this 
will prove crucial for the detection prospects once the astrophysical 
background and the angular response of the instruments are taken 
into account. 



3.3 Distribution of J factors and 030% for the cluster sample 

Most of the galaxy clusters in the MCXC are faint objects (see 
Fig. [I}. A stacking analysis is appealing if the slope of log A^ — 
log J is steeper than —1, indicating that the number of sources in- 
creases more rapidly than the brightness of those sources dimin- 
ishes. 

The logj^o (■^) ^ iogio (J) distribution is shown in the top pan- 
els of Fig. [3] We note that the double-peaked structure found is 
an indication that the MCXC is neither complete nor uniform at 
high redshift. The top left panel emphasises the importance of sub- 
structures for Qint = 0.1°: in their absence (dotted blue line), we 
have Ano-subs <x J^^'^ such that there are > 20 times more ob- 
jects each time J is decreased by a factor of ten. With substruc- 
tures (dashed red line) the prospects for stacking are improved; 
Asubs oc J^^" such that there is now a factor 100 increase in 
the number of target objects for the same factor ten J decrease. 
The lower left panel of Fig.|3]shows the logj^Q(A) — \og-^Q{agQ%) 
distribution, where Q80% is the integration angle for which 80% 
of the total J-factor is included. The quantity (ogo%) is of impor- 
tance as it corresponds to the desired PSF in order to include most 
of J in the majority of sources. This plot again emphasises the role 
of substructures. The mean for the ago% distribution moves from 
~ 0.03° (dotted blue line) to ~ 0.15° when the contribution of sub- 
structures is taken into account. This is more favourable for cuiTent 
observatories, the angular resolution of which being at best ~ 0.1°. 
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Figure 3. Top panels: the distribution of logjg A^ — log^o J(aint)- The left panel shows a comparison without (dotted blue line) and with (dashed black 



line) substractures for aint = 0.1°. The right panel shows 1 



I A'^ — logjg J(oint) for three different integration angles (all with substructures). The solid 



lines are power-law fits on the brightest J values of the histograms. Bottom left panel: the distribution of ago% C^^e integration angle containing 80% of ,7) 
without (dotted blue line) and with (dashed black line) substructures. Bottom right panel: the distribution of boost factors (for the MCXC sample) for four 
integration angles. The boost is defined to be the ratio of the total J factor (with substructures) to the J factor obtained in the hypothetical case where no 
substructures (only smooth) exist in the galaxy cluster. 



The choice for the integration angle also impacts on the 
logjQ(A'^) — logjQ(J) distribution. The top right panel shows that 
larger integration angles have an impact only for halos not fully 
encompassed, i.e., for the closest/brightest ones. Indeed, objects 
whose 080% < Clint do not have significantly more signal when 
Oint is increased. For the bigger objects the interplay between the 
different angular dependence of the smooth and substructure contri- 
butions shapes the log\o(-^) ^ logio('^) distribution. The distribu- 
tion of boost (for different integration angles) shown in the bottom 
right panel of Fig. [3] illustrates this interplay. For very small Oint 
(e.g., 0.01° dash-dotted grey line), the signal from the smooth dom- 
inates and the distribution is strongly peaked around 1 (no boost). 
As the integration angle is increased (0.5° , solid red line), the distri- 
bution is broadened, asymmetric, and reaches a maximum of ~ 20. 



3.4 Impact of varying substructure parameters 

As already underlined, several ingredients for the DM distribu- 
tions (smooth and subhalos) can affect the results above. For in- 
stance, physical processes involving baryoni c matter (such as AGN 
feedback) may produce a core distribution dMartizzi et al.ll2012f) . 
This could decrease the total J values. However, for galaxy clus- 
ters, a dominant part of the signal comes from substructures for 



Clint ^ 0.05° (as seen in bottom-right Fig. [3]l, the distribution of 
which impacts the results significantly. First, the smallest proto- 
halo mass remains unknown, and it strongly depends on the de- 
tails of the da rk matter candidate microphysics at the kinetic de- 
coupling (e.g.. Green et al.ll2005c IProfumo et al.luOOfit iBringmannI 
l2009luGondolo et al.l2012l) . Second, the subhalo spatial distribution 
is found to be less concentrated than the smooth hal o one, and con- 
sisten t with Einasto profiles in the re cent Aquarius (ISpringel et al.l 
I2OO8I) and Phoenix JGao et al.|[2012l) simulations. In the latter the 
mass distribution slope qm {dP/dM oc M^"") is also found to 
be steeper and close to 2, leading to a larger fraction of substruc- 
tures in clusters than in galaxies. When the slope is close to 2, the 
contribution to the signal of small subhalos becomes as important 
as that of larger ones, which can strongly boost the overall signal 
depending on the chosen ca — Ma (concentration-mass) relation. 
Many studies have focused on the caracterisation (e.g., mean 
and variance, environment effects) of this relation, but are limited 
by the mass resolution of cuiTently available numerical simula- 
tions. The state-of-the-art studies on galaxy clusters apply down 



to a minimum ha l o mass ~ 10 Mp^ dWechsler et al 
Zhao et all I2OO3I |2009|: iNeto et aUHool iGao et alT 



ga 



apply 

20021 



2011, 



Maccio et a l.'2008VGi ocoli et al.l201(]| ; lMunoz-Cuartas et al.l2011 



2ooe 



200s 



Klvpin et al. 2011 ). Extrapolations to the smallest subhalo mass are 



provided in a very few studies only. For instance, almost all anal- 
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yses of the DM annihilation signal are based on tw o different pa- 
ramete risations ( Bullock et al. 2001 ; Eke et al. 2001 ). lGiocoli et alj 
( 120120 recently provided a n ew parameterisation of ca — Ma : at 
z = 0,it is consistent with the lBullock et alj ( l200in parametrisation 
in the < 10^'' M0 range, but its redshift dependence is different. 

Given the variety of results found in the literature and the un- 
certainties on some parameters, we only select a few configurations 
below. 

• dV / dVphocnix' uses spatial distribution and scale radius as 
provided by the Phoenix project dOao et al.ll2012h . instead of fol- 
lowing that of the smooth profile (all other parameters as in Sec- 
tion|2J2j; 

• dP/dVphocnix and om ~ 1.98: uses a ste e per sl ope for the 

lOao et al.l I2OI2I) instead of 



mass distribution as found in Phoenix 



1.9 (note that au = 1.94 in lSpringel et"ai]|2008h ; 

• dP/dVphocnix, Qm = 1. 98, and / = 0. 3: uses a DM mass 
fraction as found in Phoenix dOao et al.ll2012h instead of 0.1 (as 
found in lSpringel eTaLll2008h . 

• dP/dVphoenix, am = 1-98, / = 0.3, and (ca — 
AfA) ENS0i or (ca — AfA) Gi2: uses a mass-con centration relation 
from lEke et al.l (1200 ih and iGiocoli et al.l ( l2012h . instead of using 
lBullocketal.l l l200lh . 



The impact is shown in Fig. |4] for the logj^Q(Af) — logj^p(J) (top 
panel) and the logj^Q(A'') — log j^q (Boost) (bottom panel) distribu- 
tions for Qfint = 0.1°. Taking an Einasto profile for the spatial dis- 
tribution of subhalos has only a minor effect (solid thin- vs dashed- 
black line). A major impact is that of the value of the parameter 
um (solid thin- vs dotted-black line). It increases the boost (and 
thus the signal) by about one order of magnitude. This increase 
can be larger if a smaller minimal mass for the subhalos is cho- 
sen, but it can also be decreased by a factor of ten if the minimal 
mass allowed is W^Mq. Although the Phoenix and Aquarius sim- 
ulations tend to prefer values close to 2, the result is intrinsically 
limited by the mass resolution of the simulation, and this slope can 
still be smaller at lower masses. Moreover, |Elahi et al. (2009) ar- 
gue that this slope could be overestimated, even in the simulation 
mass range. Another obvious effect is from the mass fraction / in 
substructures (thin vs thick solid black line). Finally, the two red 
curves show the impact of the ca — Ma p arametrisation (thin and 
thick red lines) compared to using iBullock et al. (20011) para metri- 
sation (thick solid black line). Ch oosing Giocoli et alj ( 120120 gives 



slightly less signal, whereas using lEke et al.l ( 1200 ih washes out the 
boost completely. To conclude, we see that the main source of un- 
certainties correspond to the slope of the mass function, the min- 
imal mass of the subhalos, and the concentration. Unfortunately, 
these are also the least-constrained parameters from the available 
numerical simulations. 

Our reference configuration gives a conservative estimate of 
the signal expected from galaxy clusters. We underline that the 
shape of logj^g A'^ — log^Q J(aint) is only weakly impacted by 
choosing other configuration (and the ordering of the best tar- 
gets — not shown — also remains mostly unaffected). Therefore, the 
conclusions that will be reached below for the stacking analysis us- 
ing the reference configuration hold regardless of this choice. How- 
ever, the consequences of a larger boost would be the following: i) 
a larger signal and a better (av) limit set from non-detection, ii) 
an increase of the 80%-containing angle, iii) an enhanced contrast 
with respect to the Galactic background resulting in an increased 
extragalactic to Galactic signal ratio. 
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Figure 4. The distribution of logjg A'^ — logio J{ctlnt) (top panel) and 
boost factors (bottom panel) for different substructure configurations (see 
text for details). 



4 HALO STACKING AND RESULTS 

There are two primary considerations for the MCXC stacking anal- 
ysis: how to order the sources, and how many sources to stack. This 
is discussed for different situations before moving to the detection 
prospects for the stacking strategy. 



4.1 Strategy for a 'perfect' instrument 

Signal-limited regime The top panel in Fig. [5] shows the cumula- 
tive distribution of J for integration angles of 0.05° (solid-green 
stars), 0.1° (open-black squares), and 0.5° (solid-red circles) as 
well as ago% (open-blue circles). The numbers denote the num- 
ber of MCXC contributing to the cumulative in a given J bin. The 
MCXC sources are naturally ordered by J in this plot. Sources 
within 20° of the Galactic centre are excluded. The cumulative 
Joai is also shown (dashed lines). As mentioned in Section 13.21 
the contrast ( Jtargct/Jcai) is related to the detectability of an ob- 
ject if we are only limited by the amount of signal available. In 
such a regime a stacking analysis remains valid as long as we 
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central part of halos, it is subject to 7-ray and cosmic-ray contam- 
ination from astropliysical sources. Among these are the power- 
ful AGN (hosting a super-massive black hole) o ften found at the 
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Figure 5. Top: The cumulative J (i.e. Y2i Ji f™ ^'1 * f°'' which Ji > J). 
The signal and associated Galactic DM background are represented by an 
arrow and a line respectively. The open-blue circles correspond to an in- 
tegration angle for which 80% of the total J of a galaxy cluster is in- 
cluded. Three integration angles are shown: 0.5° in sohd-red circles, 0.1° 
in open-black squares, and 0.05° in solid-green stars. Clusters closer than 
0;.^ = 20° from the GC are discarded. Bottom: J]i Ji{> ^)/«intV^ 
(proportional to the cumulative signal-to-noise for a fixed integration angle) 
as a function of J for the same integration angles. 



add sources with a contrast larger than one. The boxed number 
in italics indicate at what point this occurs: for an integration an- 
gle of 0.5° the optimum number of objects to stack in this regime 
is 21. The wealth of sources in the MCXC becomes more useful 
for smaller integration angles, with an optimum of 1224 objects 
at 0.05° . For the latter, the contrast never falls below one, but be- 
yond 1224 objects, the total J does not significantly increase. For 
aso% only 10 sources can be stacked before the signal is domi- 
nated by Jcai- The total J (with a contrast> 1) available in these 
scenarios is Yl ^numbered « 4 X 10^^ Mq kpc~^ for Qint — 0.5°, 
8 X 10^^ Mq kpc"^ for aint = 0.1°, and 5 x 10^^ Mq kpc"^ 
for Qint ~ 0.05°. The maximal value that can be achieved is 
J2 "^numbered ~ 2 x 10^^ Mq if Oint = 080%. i-S- ten times the 
result that can be achieved at fixed integration angle. 

Background-limited regime: all-sky vs pointed instruments It 

is not just the Galactic DM background that is important in the 
selection of target objects, but also the astrophysical 7-ray back- 
ground. As the DM annihilation signal is prominent at the very 



2007b. or intra-cluster 



cluster centre (e.g., JMcNamara & Nulsen 

shock-driven particle acceleration (e.g., J Markevitch & VikhlininI 
I2OO7I : lEnBlin et al.l I2OIII : IPinzke et al]|201ll) . This astrophvsical 
background will increase with the square of the integration angle. 
The signal-to-noise ratio for a source is therefore proportional to 
J/\/{a ). The cumulative signal-to-noise ratio for an all-sky in- 
strument (in which all objects are observed for the total observa- 
tion time) is therefore proportional to J2i Ji{> J)/ Yl ^i- Fo'' ^ 
fixed integration angle, this is J^i >^i ( > ^) /Qint ^/N . For an instru- 
ment that relies on pointed observations, the amount of time spent 
on each source is the total observing time available divided by the 
number of sources that must be observed. Therefore, the signal-to- 
noise ratio is proportional to J]]^ Ji (> J)/Y1 ctfVN. In that case, 
the best strategy appears to focus on a single bright object. As the 
total available observation time is fixed, time spent observing addi- 
tional sources reduces the time spent observing the brightest target 
(see Section |431 l. 

The lower panel of Fig. |5] shows X]i Ji{> J)/aintVN as a 
function of J, again for integration angles of 0.05° (solid-green 
stars), 0.1° (open-black squares), and 0.5° (solid-red circles). The 
peak in these 'signal-to-noise' curves indicates the optimum num- 
ber of sources to stack in the background-limited regime, and are 
highlighted as 1224, 713 and 21 for 0.05°, 0.1° and 0.5° respec- 
tively. In this plot, sources are ordered by increasing J-values, and 
therefore only 'signal-to-noise' curves can be included for fixed in- 
tegration angles. For variable integration angles, such as a8o%. the 
signal-to-noise ratio of each source in the catalogue will depend on 
the integration angle as well as J, and therefore the stack must be 
ordered by J/ago% • In this case the optimum number of sources 
is close to the full stack size, though we will see in the following 
section that these optimum values change drastically when the an- 
gular response of the instrument is considered. Examining the list 
in detail, it is apparent that when ordering by J/ago% rather than J, 
only a few sources high-up the list swap places. The sources falling 
somewhere in the 'top' 20-30 remain consistent. 

The conclusions drawn from Fig. [5] are only valid for a in- 
strument with a perfect angular response. In reality, the angular 
response of an instrument — typically characterised by the point 
spread function (PSF) which we take here to mean the 68% con- 
tainment radius — must be combined with the integration angle in 
quadrature before considering the amount of background contami- 
nation in an observation. In deciding which integration angle to use, 
we consider that, for a small fixed angle, the cumulative J is re- 
duced since some signal from angularly-large sources is neglected. 
For a large fixed angle (e.g., 0.1°), the cumulative J increases 
slowly, implying that angularly-large sources are also bright, and 
located near the top of the list. Further down the list, where sources 
are angularly small, large amounts of galactic contamination and 
astrophysical background are included unnecessarily. Therefore a 
different integration angle for each source, such as ago%, may be 
optimum, and is used in the remainder of the analysis. 



4.2 Strategy for a 'real' (PSF-limited) instrument 

The upper panel of Fig.|6]shows the cumulative signal-to-noise ra- 
tio as a function of the number of sources stacked for different val- 
ues of the PSF. As the PSF worsens from 0.01° to 3°, the relative 
signal-to-noise ratio drops and the peak position shifts towards a 
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Figure 6. Top panel: The signal-to-noise ratio as a function of the number 
of sources to stack for different values of an instrument point-spread func- 
tion (PSF). Bottom panel: Optimum number of sources as a function of the 
PSF for Oint set to a8o%- For a fixed integration angle ajnt = 0.1°, this 
number is constant with PSF (dashed line). 
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Figure 7. Top panel: the 5o" sensitivity of Fermi-LAT (5 years exposure) 
(solid curves) and CTA (1000 hours exposure per source) (dashed curves) to 
the three brightest MCXC sources in the context of this work: Virgo (black), 
A 426 (red) and A 3526 (blue) when considering an integration angle of 
agO% . Bottom panel: as above, for stack sizes of the optimum number of 
sources for a 0.1° (1200) (blue), 0.5° (90) (green) and 1° (17) (red) PSF 
obtained from Fig. |6] Virgo alone is again shown in black. For CTA the 
1000 hour exposure is divided equally over the number of sources in the 
stack. 



smaller stack size. The peak position indicates the optimum num- 
ber of sources to stack, and is shown in the lower panel of Fig. |6] 
as a function of PSF for an all-sky instmment. For a fixed inte- 
gration angle of 0.1° (dashed line), the optimal number is constant 
with the PSF. When 030% is considered, the optimal number of 
sources drops as the PSF of the instrument increases. For a PSF 
of O.r, 1200 sources should be stacked. For a PSF of 0.5°, 90 
sources should be stacked and for a PSF of 1 ° , 17 sources should be 
stacked. When the PSF increases above ~2°, stacking is no longer 
a valid approach, and only the brightest source should be consid- 
ered. It is not only the number of sources that should be stacked that 
changes with PSF, but also the order of those sources. Independent 
of the PSF the top two sources are Virgo and then A 426. At a PSF 
of >1° the third brightest source is NGC4636. However, below a 
PSF of 1°, A 3526 moves into third place. The top ten sources al- 
ways contains Coma, but Fornax falls out of the top ten when the 
PSF drops below ~0.1°. 



4.3 Detection Prospects 

In this section, we assess the DM detection prospects for the 
stacking of sources from the MCXC for the Fermi-LAT all-sky 
7-ray satellite, and the envisaged array of Imaging Atmospheric 
Cherenkov Technique CTA (Cherenkov Telescope Array). Whilst 
the design of CTA is still evolving, performance curves for several 
configurations have been release d. Here, we use the so-cal led array 
layout 'E', which is described in lCTA Consortium! (120101) . For the 
Fermi-LAT, the 1-year poi nt-source performan ce curves for a high- 
latitude source are used dRando et al.l l2009h . The diffuse galac- 
tic and extra-galactic background models given by the template 
files gal_2yearp7v6_v0 . fits and iso_p7v6source . txt 
respectively, which are available from the Fermi-LAT data server, 
are used to obtain the background within the integration angle for 
each so urce position on the sky. A toy likelihood-based model, as 
used in lCharbonnier et al.l ( 1201 If) , is used to obtain the sensitivity 
of these instruments to the DM galaxy cluster signal. 

The top panel of Fig.|7]shows the sensitivity of Fermi-LAT and 
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CTA to the three sources from the MCXC that result in the high- 
est J and J/ajio% ('signal-to-noise') in the context of this work 
for PSFs smaller than 1°: Virgo, A 426, and A 3526 when consid- 
ering an integration angle of a80%- All curves represent a 5-sigma 
significance. The Fermi-LAT curves are computed for 5 years ex- 
posure, whilst the CTA curves assume 1000 hours observation of 
each source. Whilst it is unrealistic to expect pointed observations 
on all three of these sources for this duration, an equal exposure is 
useful in comparing the potential targets. Virgo dominates the sen- 
sitivity for both detectors. Individual curves were produced for the 
ten brightest sources, and it was found that the top three sources 
shown here provide the best individual sensitivity. 

A spectrum of photon energies is associated with each DM 
mass. Most sensitivity is contributed by the photon energy range 
close to the peak in E^dN/dE, which lies one order of magnitude 
below the DM mass for our assumed annihilation spectrum. Very 
low energy photons (several orders of magnitude below the DM 
mass) contribute little to the sensitivity due to the relatively hard 
signal spectrum and overwhelming background. In our analysis we 
exclude photons with energies less than 1/200 the DM mass (pro- 
viding this cut lies below 10 GeV). We consider this to be a realistic 
approach in practice to avoid source confusion problems due to the 
very poor PSF of Fermi-LAT close to threshold. At 100 MeV for 
example, the Fermi-LAT PSF is some 6° in radius, a region that is 
likely to include several additional Fermi sources. 

The lower panel of Fig.|7]shows the sensitivity of Fermi-LAT 
and CTA when stacking the optimum number of sources deter- 
mined from the lower panel of Fig. |6]f or PSFs of: 0.1° (1200), 0.5° 
(90) and 1° (17). The brightest source, Virgo, is shown individu- 
ally. Again, the Fermi exposure is taken as 5 years. As Fermi is an 
all-sky instrument, each source in the stack receives this exposure 
regardless of the stack size. 

At DM masses below ~100 GeV, the majority of photons are 
collected in an energy range where the Fermi-LAT PSF is worse 
than a degree. Here, the analysis falls into the background-limited 
regime. Therefore stacking does not help, and just adds back- 
ground, making the sensitivity worse than Virgo alone, for exam- 



ple: ~3.5xlO"^° cm'' s" 



to 



-2.5x10 ^^ cm^ s S-espectively 



at ~2 GeV. Searching for WIMP masses above ~100 GeV, pho- 
tons begin to be included that are seen by Fermi-LAT with a better 
PSF. In this mass regime, the amount of signal collected becomes 
important, and the stacking helps. The analysis eventually becomes 
signal-limited, and stacking improves the sensitivity by a factor of 
up to 1.7, from ~3xlO"^^ cm^ s"^ to ~1.8xI0"^^ cm^ s"^ at 
~I TeV. This is roughly equivalent to the improvement in signal- 
to-noise ratio shown in the upper panel of Fig. |7]for a PSF repre- 
sentative of the energy range in question. For example, at a mass 
of 2 TeV, photons are included down to 10 GeV, corresponding to 
a PSF always better than 0.25° . Even at masses where an improve- 
ment with stacking is found, beyond a stack size of 17 sources the 
improvement is negligible. This is simply because the instrument 
PSF varies with energy and therefore taking the optimum number 
of sources for a fixed PSF is only an approximation. 

In the case of CTA, we assume that a total exposure of 1000 
hours is available, and since CTA requires pointed observations, 
this is reduced to ~60 hours per source when 17 objects are 
stacked, ~12 hours per source when 90 objects are stacked, and 
~0.8 hours per source when 1200 objects are stacked. This ef- 
fect dominates any gain in sensitivity due to stacking, and confirms 
the finding of the previous section that for an instrument requiring 
pointed observations, only the brightest source should be targeted. 



Note that systematic effects are not included, and will limit the ac- 
curacy of a 1000 hour observation. 



5 DISCUSSION 

A stacking analysis of galaxy clusters may provide better limits for 
indirect detection of DM than the analysis of any single object, at 
least for all-sky instruments. However, this improvement is likely to 
be modest for the case of annihilating dark ma tter Stacking is mor e 
promising in the case of decaying dark matter dCombet et al.l2012h . 
For instruments requiring pointed observations such as CTA, ob- 
serving the most promising source until the observation is system- 
atics limited and then moving to additional sources is a reasonable 
strategy. Such an approach also mitigates against the uncertainty in 
the properties of individual halos. 

Limits placed on the velocity-averaged cross section depend 
on the determination of J not only for studies relying on known 
detector sensitivities (such as this work), but also for works mak- 
ing use of real data, e .g. the Fornax observation by H.E.S.S. 
JAbramowski et al.ll2012h . We checked that given the same J for a 
given source, we obtain a very similar sensitivity to that estimated 
in previous studies (see Appendix |B]l. In our analysis, Virgo has 
the highest astrophysical factor (J) and best signal to noise ratio, 
followed by A 426. Several authors have suggested (based on clus- 
ter properties given by the HIFLUGCS catalogue) that Fornax is 
the most promising galaxy cluster for DM annihilation. However, 
as discussed above, the MCXC provides homogenised values for 
A/500 based on a more accurate gas density prescription that typ- 
ically results in lower J for the brightest clusters (but note that 
the re is no systematic tre nd when all galaxy clusters are compared, 
see lPiffaretti et al.ll20I lb . The differences between these two cata- 
logues are large enough to significantly change the conclusions of 
studies on the sensitivity of current and future instruments to DM 
annihilation, for example the detectability (or not) of DM with the 
annihilation cross-section expected for a thermal relic in this class 
of objects. In that respect, the ranking we provide from the MCXC 
catalogue should be robust, although the J values calculated in this 
paper may still change depending on the level of dumpiness, exact 
mass-concentration relation, etc. 

For all-sky instruments and in particular for Fermi-LAT, the 
improvement in sensitivity obtained by stacking is at best a factor 
1.7: MCXC sources with the 1200 largest values of J or J/aso% 
should be included to obtain this improvement. Additional sources 
do not improve the sensitivity, as further background is integrated 
without significant additional signal. This implies that the benefits 
of stacking are limited by the PSF of the available all-sky 7-ray 
instruments. Indeed, the PSF of Fermi-LAT at low energy is sev- 
eral degrees, while the majority of MCXC targets are distant and 
hence subtend small angles, with a typical Q80% of ~ 0.15° (when 
substructures are considered): an all-sky instrument with a PSF ap- 
proaching UgQX at all energies would benefit from the stacking of 
all sources in the MCXC. In this case, sensitivity would then be 
limited only by the available signal, and an extended catalogue — as 
shoul d be provided in a f ew years from now by the eROSITA mis- 
sion iPredehl et al.ll20I ih — including even fainter objects would 
be needed to reach a cumulative J ~ 10^^ — 10^^. 

A stack of the top 1200 objects excluding Virgo results in 
a sensitivity only ~15% worse than the same stack size includ- 
ing Virgo. In this case, the improvement in sensitivity between 
the brightest source alone (A 426) and the stack of 1200 objects 
is nearly a factor of 3 above masses of 100 GeV. The advantage is 
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that the large number of clusters stacked is expected to wash out in- 
dividual uncertainties on the halo properties (e.g., the dispersion of 
mass-concentration relationship). One viable strategy might there- 
fore be to use Virgo as an independent confirmation of the signal es- 
tablished through the stack ing of other clusters. Virgo contains th e 
known 7-ray emitter M 87 dSeilicke et al.ll2004lAbdo et al.ll2009l) . 
The ago% of Virgo is ~0.3° for a smooth halo, comparable to the 
Fermi-LAT PSF at the highest photon energies, but ~3° when sub- 
structure is considered. Disentangling the point-like emission from 
M 87 from an y extended DM e mission may therefore be possible. 
Very recently, iHan et al.l l|2012i) have claimed evidence at the ~4(t 
level for diffuse DM-like emission from Virgo: they use photon en- 
ergies detected by Fermi-LAT above 100 MeV and a full likelihood 
fit to a template vs. a point source. Further Fermi-LAT observa- 
tions, and deeper investigation of possible astrophysical origins for 
the apparent extended emission, are required to confirm or refute 
this intriguing result. 

The great advantage of an all-sky instrument such as the 
Fermi-LAT is the simultaneous observation of all sources. Anal- 
ysis of the potential DM signal from galaxy clusters can therefore 
be performed for different numbers of stacked objects with differ- 
ent orderings simultaneously. In the event of any detection from a 
stacked analysis, a re-analysis on a different, more numerous, set of 
objects may help to confirm the result. CTA only becomes compet- 
itive with Fermi for DM masses above ~ 1 TeV. However, at these 
energies CTA will have an angular resolution approaching 0.02° 
and may therefore help in isolating point-like sources from clusters 
(Virgo may not be the only galaxy cluster with a 7-ray emitting 
source embedded within), to aid in the choice of sources to stack 
for a Fermi analysis, or in a hopeful case to rule out a point-like 
emitter as the source of Fermi detection. CTA may also be critical 
to measure the cut-off in the DM annihilation spectrum for heavy 
dark matter, and hence measure the DM mass and establish the uni- 
versality of the annihilation spectrum. 

Data analysis can be optimised by adapting the integration 
region for each cluster, as we have shown with the example of 
0^80% • We provide the necessary ingredients to refine the anal- 
ysis presented here in Appendix |A] From the dark matter mod- 
elling side, a systematic study remains to be done to take into ac- 
count various DM profiles, substructure characteristics, the mass- 
concentration dispersion, etc. This will be carried out in a future 
work. We reiterate here that our limit on (av) could be changed by 
taking other configuration of the substructure distribution. In the 
most favourable case, it would allow to reach the benchmark value 
(av) ~ 3 lO"^*" cm^ s~^ coming from cosmological constraints. 
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APPENDIX A: THE RELATIONSHIP BETWEEN J AND 

aiNT 

There exists a simple para metrisation to calcul ate J(aint) for any 
ctint, given the DM profile JCombet et al.l20I2r) . Indeed, we can as- 
sume that all galaxy clusters share the same DM profile. Given the 
mass range span by the MCXC, we can approximate at first order 



their concentration parameter to be the same. For a NFW profile , 
c{M) = i?vir/r-s and we take c(1O"M0) ~ 5 JDuffv et alJJOOsh . 
Defining 



tan 



" (7) ' "" 



; tan 



1 / 5rs 



Qint J 

X = , and Xn 

Cts 



(Al) 

(A2) 



there is a universal depen dence of the fractio n of the smooth and 
substructure contributions (IMaurin et al.ll2012h . 



Tjix) = 



,J{x ■ as) 
7 ' 



which we parametrise to be (valid only for a NFW) 

(3x°'>^ ifa;<10-2, 
^smooth (a;)=< 1 if X ^ 5, 

[-0.086+0.17 lii(a:)-0.092 In^ (2!) + 0.011 ln^{2^)] . 



and 



-^subs(a;)= 



if 3; ^ 5, 

[-1.17+1.061n(a:)-0.171n^{2!)-0.0151n^(a;)l 



(A3) 



(A4) 



The 'signal' J can then be calculated for any integration angle. 



using 

Jtot(aint) — Jsmooth(0-1°) X 
+ Jsubs(0.1°) X 



-^sm (Qjnt/Qs) 

j;m(0.1°/a.) 

-^subs (aint/cMs) 



(A5) 



j;ubs (0.1° /a,) ' 

Hence, as shown in lMaurin et alj l l2012n . for DM annihilation, one 
needs three quantities (available for all clusters in the Supplemen- 
tary Material — ASCII file — submitted with the paper, short sample 
in appendix O. i-e-> cts, Jsmooth(0.1°) and Jsubs(0.1°). 

This parametrisation describing the fraction of the signal in a 
given angular region is valid down to J^j = 10^'*. 



APPENDIX B: A COMPARISON OF THE VALUES OF J 
OBTAINED HERE TO OTHER WORK 

DM annihilation in gal axy clusters ha s been studied in sev- 



eral papers including (Jeltema et al.ll2009 : 
Sanchez-Conde et al.l 2011: Pinzke et al 



Ackermann et al.l 



201 ll: iGaoetalJ 



201c ; 



2012 : 



Huang et alj|2012l : lAndo & Nagaill2012l : Man et alj|2012l) . Below, 
in Table IbT] we provide a comparison with some of these studies, 
whenever the J factor was available. 

The calculations of the pr e sent work are consis tent with those 
oflSanchez-Conde et all ( 1201 ih . lPinzke et alj ( I20TTI) and lOao et al] 
I2OI2I Our results for the boost values are also in agreement. Our 
J values are also broadly consistent tho ugh systematically lower 



all y lower 
ai1l l20I(]|) 



(res p. systematically lar ger) than those of lAckermann et 

and lHuang et al.l j2012r) without (resp. with) the substructure con- 
tribution. In any case, the uncertainties quoted in these two papers 
(third line in the table) is probably underestimated. Note that all the 
three studies rely on the HIFL UGCS catalogue based on ROSAT 
and ASCA X-ray observations (IChen et al.ll2007h . The main differ- 
ence is for Fornax, which is a factor of ten l arger (though the dif- 
feren ce is less significant if we compare with lSanchez-Conde et al.l 
I2OIII results). This is due to the lower mass we infer for this clus- 
ter from the MCXC il/soo. and i?5oo val ues, which are based on a 
better modelling of the gas in the cluster ( IPiffaretti et al.l20I ih . 
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Table CI. J values of the five first objects of the MCXC catalog. The full table is available in the onhne version. 



NAME 



Indice MCXC 



I 
[deg] 



b 
[deg] 



d 
[kpc] 



013 
[deg] 



Jsm(O.ldeg) 

[Af|,kpc-5] 



Jsub(O.ldeg) 
4^ 



[Af2kpc-5] 



"80% -^80(^80%) 



[deg] 



[M|,kpc-5] 



UGC 12890 


1 


RXCJOOOO.4.0237 


2 


RXCJOOOl.6.1540 


3 


A2692 


4 


A2697 


5 



101.78 -52.48 1.62e-l-05 8.36e-02 

94.27 -62.62 1.55e-i-05 6.23e-02 

75.13 -73.73 4.61e-l-05 4.01e-02 

104.31 -49.00 6.89e-i-05 3.19e-02 

92.17 -66.03 7.63e-l-05 3.77e-02 



9.69e-i-08 


8.05e-l-09 


2.73e-01 


1.94e-HlO 


5.31e-l-08 


5.32e-l-09 


2.13e-01 


l.Ole-l-10 


2.82e-l-08 


3.35e-i-09 


1.27e-01 


4.30e+09 


2.16e-l-08 


2.41e-l-09 


l.OOe-01 


2.67e+09 


3.30e+08 


3.19e4-09 


1.14e-01 


3.88e+09 



Table Bl. Comparison with J values from other works for a-^ 
and 1° respectively. 



0.1° 





log 


( '' 


(1°) 


■) 


logio 


( J(0.i°) ^ 




LU \GeV- 


'- cm-5 


VGcV^ cm -5 ) 


Ref. 


[1]* 


[2]§ 


This work 


[3]^ 


This work 


Error 


SO.l 


<0.2 


(wo/w subs) 


- 


(no subs) 


Fornax 


17.8 


17.9 


16.9 


18.8 


17.0 


16.7 


Coma 


17.2 


17.1 


16.9 


18.4 


16.8 


16.7 


A1367 


- 


17.1 


16.7 


18.3 


- 


16.5 


A1060 


- 


17.3 


16.8 


18.3 


- 


16.7 


AWM7 


17.1 


17.2 


16.8 


18.2 


- 


16.6 


NGC4636 


17.6 


17.5 


17.2 


18.2 


- 


16.9 


NGC5813 


- 


17.3 


17.1 


18.1 


16.4 


16.8 


A3526* 


17.4 


- 


17.1 


18.1 


- 


16.9 


A426t 


- 


- 


17.2 


18.1 


16.9 


17.0 


Ophiuchus 


- 


- 


16.8 


18.1 


16.8 


16.7 


Virgo 


- 


- 


17.9 


18.0 


17.5 


17.5 


NGC5846 


- 


- 


16.7 


17.9 


16.5 


16.5 



[kckermann et alJi201(]|)jHuang et ai]i2012l) . 1Sinchez-Conde et al.N201lh 
*Centaurus, tperseus 

APPENDIX C: J TABLE OF THE MCXC CATALOG 

In Table |CT] the J values of the five first objects of the IVICXC 
catalog are given. 

We provide in the online version the full table for the 1743 IVICXC 
catalog objects. 

Table IC I [ column description : 
NAJVIE: Cluster name (from IVICXC) 

Indic e MCXC: Indice of the cluster row in the IVICXC catalogue 
(see ( |PiffarettietaLll20Ilh ). 
I: Galactic longitude in degree (from MCXC). 
b: Galactic latitude in degree (from MCXC). 
d: Angular diameter distance in kpc (from MCXC). 
Os : arctan (rs/d) in degree (rs is the scale radius of the cluster). 
Jsm(O.Ideg): Astrophysical annihilation term (from smooth halo) 
in [Mgkpc"^] for aint = 0.1 deg. 

Jsub(O.ldeg): Astrophysical annihilation term (from substructures) 
in [M|kpc"''] for aint = 0.1 deg. 
080% ■ Angle in degree containing 80% of the total J. 
J80: Value of annihilation signal for aint = a80%- 
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